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For the filtering of peaks in periodic signals, we specify polynomial filters 
that are optimally localized in space. The space localization of functions / 
having an expansion in terms of orthogonal polynomials is thereby measured 
by a generalized mean value s(f). Solving an optimization problem includ- 
ing the functional we determine those polynomials out of a polynomial 
space that are optimally localized. We give explicit formulas for these op- 
timally space localized polynomials and determine in the case of the Jacobi 
polynomials the relation of the functional e(f) to the position variance of a 
well-known uncertainty principle. Further, we will consider the Hermite poly- 
nomials as an example on how to get optimally space localized polynomials 
in a non-compact setting. Finally, we investigate how the obtained optimal 
polynomials can be applied as filters in signal processing. 
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1 Introduction 

For the detection of peaks in mass spectrometry data, window functions are often used 
to preprocess the incoming signals, for instance, to perform a baseline correction or to 
filter out disturbing higher frequencies (cf. |22j . |29j). In the following, we consider 2n- 
periodic signals / filtered by a convolution with a window function h, i.e. the filtered 
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signal Fhf is given by 



F h f{t) := (f *h)(t) = — f(s)h(t - s)ds 



i r 



(i) 



For peak detection purposes, two properties of the window function h are important. 
First of all, the function h should be localized in frequency in order to filter out the high 
frequencies. On the other hand, the window h should also be well localized in space 
such that the convolution operator Fh is still able to determine the peaks of the signal 
/. However, the uncertainty principle states that it is impossible that the function h is 
well-localized both in space and frequency (see, for instance, [3], [7], [IH])- Therefore, 
in search for an optimal filter h for peak detection, one has always a trade off between 
denoising the signal / and determining the position of the peaks of /. 

The main objective of this article is to investigate the space localization of polynomial 
filters. If the window function h is a trigonometric polynomial of degree n, the frequency 
domain of h is bounded. Therefore, a polynomial filter h in a peak detection process 
automatically performs a low-pass filtering of the signal /. It remains to analyze the 
space localization of h. In this regard, one has to specify how space localization of a 
function / is measured. 

Beside the trigonometric setting, we will consider general systems of orthogonal polyno- 
mials in this article. For functions / having an expansion in orthogonal polynomials, we 
will use a functional s(f) to measure the space localization of /. We will study an opti- 
mization problem including this functional s(f) that allows us to construct polynomials 
and band-limited functions that are optimally localized in the space domain. In partic- 
ular, for Jacobi polynomials, it will turn out that the functional e(/) is related to the 
position variance varg(/) of a well-known uncertainty principle and that the optimally 
space localized polynomials minimize also the term for the position variance vais(/). 

By now, research in this direction has been done mainly for trigonometric polynomials 
on the unit circle and spherical harmonics on the unit sphere E> d . In [24j, Rauhut used 
the position variance of the Breitenberger uncertainty principle to construct optimally 
localized polynomials on the unit circle. On the unit sphere, the works of Mhaskar et al. 
|20| and Lain Fernandez |15j led to optimally space localized polynomials and polynomial 
wavelets on the unit sphere S d . One aim of this article is to extend these results to 
orthogonal expansions on the interval [—1,1]. These more general results are then used to 
construct new window functions for peak filtering purposes in the trigonometric setting. 

In the following, we will mainly consider the Hilbert space L 2 ([-l,l],w) with the inner 
product 



J-l 

where the weight function w is a nonnegative continuous function on [—1, 1]. Then, for 
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f € L 2 ([— 1, 1],u>)j we define the functional e(f) by 

-l 

|2„ 



8(f) := J ^x\f{x)\ z w{x)dx. (2) 

Exactly this generalized mean value e(f) is the starting point for our investigations. In 
the following first section, we will see that the mean value s(f) determines a measure for 
the localization of the function / on the boundary values of the interval [—1,1]. Then, we 
will study the functional e(/) for polynomial subspaces of L 2 ([— 1, 1], w). In particular, 



in Theorem 2.8 and in Corollary 3.2, we will give representations of those polynomials 
V n that maximize e(/), i.e. those polynomials that are optimally localized at the right 
hand boundary of the interval [—1,1] with respect to e(f). 



it 



In the third section, we will emphasize on the Jacobi polynomials. In Theorem 4.3 
will be shown that the position variance vars(/) of an uncertainty principle for Jacobi 
expansions is minimized by the polynomials V n . 

In the fourth section, we will consider the Hermite polynomials on the real line and 
determine optimally space localized polynomials in this setting. The obtained results 
are an example of how it is possible to generalize the theory of Section [2] and [3] to a 
non-compact setting. 

Finally, in the last section we will turn back to the peak filtering application mentioned 
at the beginning. We will investigate how the optimally space localized polynomials of 
this article are related to well-known polynomial window functions and how they can be 
applied as filters for the peak detection of a signal. 



2 Optimally space localized polynomials 



We start out by introducing particular polynomial subspaces of L 2 ([— 1, 1], w). Therefore, 
we denote by {pz}^ the family of polynomials that are orthonormal on [—1, 1] with 
respect to the inner product (-,-) w . Further, we assume that the polynomials pi of degree 
I are normalized such that the coefficient of x is positive. Then, the family {pz}^ 
defines a complete orthonormal set in the Hilbert space L 2 ([— 1, 1], w) (cf. |28[ Section 
2.2]). 

Definition 2.1. As subspaces of the Hilbert space L 2 ([— 1, 1], w), we consider the fol- 
lowing three polynomial spaces: 

(1) The space spanned by the polynomials pi, < I < n: 

II n := jp : P(x)=J2ci P i(x), co,...,c„ec|. (3) 
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(2) The space spanned by the polynomials pi, m < I < n: 

II™ := i P : P(x) = C ^), Cm, ■ ■ ■ , Cn G C 1 . (4) 



l=m 



m—1 



(3) The space spanned by a polynomial 7£(x) = p m {x) + eipi{x) of degree m and 
the polynomials p/, m + 1 < Z < n: 



1=0 



U%:=lP: P( x ) = c m K(x) + ^ qpiO*), c m , . . . ,c n £ C \ . (5) 

L l=m+l ) 

Further, we define the unit spheres of the spaces U n , II™ and II? as 





= {p s n n 


||P| 




= 1} 




= {p g n™ 


: ||P 


1 10 


= 1} 






: ||P 


I w 


= 1} 


Remark 2.2. Clearly, II™ C U n 


and C 


n n . 


In 


the 



•ature, the spaces H™ are 
sometimes called wavelet spaces and considered in a more general theory on polynomial 
wavelets and polynomial frames, see |21| and the references therein. For special choices of 
1Z, the polynomials in the spaces il^ play an important role in the theory of polynomial 
approximation. In particular, if polynomial reproduction is requested, a common choice 
for the polynomial 1Z is the Christoffel-Darboux kernel of degree m (see [6], |19]). The 
standardization e m = 1 for the highest expansion coefficient of the polynomial 1Z causes 
no loss of generality and is a useful convention for the upcoming calculations. 

The first goal of this section is to study the localization of the polynomials in the spaces 
Il n , II™ and at the right hand boundary of the interval [—1,1] and to determine 
those polynomials that are in some sense best localized. As an analyzing tool for the 
localization of a function / G L 2 ([— 1, 1]) at the point x = 1, we consider the mean value 
as defined in (g. If ||/|| tt = 1, then -1 < e(/) < 1, and the more the mass of 
the L 2 -density / is concentrated at the boundary point x = 1, the closer the value e(f) 
gets to 1. Therefore, the value e(/) can be interpreted as a measure on how well the 
function / is localized at the right hand boundary of the interval [—1,1]. We say that / 
is localized at x = 1 if the value e(/) approaches 1. 

Now, our aim is to find those elements of the polynomial spaces II n , n™ and that 
are optimally localized at the boundary point x = 1. In particular, we want to solve the 
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= arg 


; max e(P), 

PeS„ 


<pm 
' n 


= arg 


; max e(P) 

Pes™ 


V n 

1 n 


= arg 


; max e(P) 

P6§? 



following optimization problems: 

(6) 
(7) 
(8) 

Remark 2.3. Since the linear spaces II n , II™ and 11^ are finite-dimensional, the unit 
spheres S n , §™ and §^ are compact subsets and the functional e is bounded and contin- 
uous on the respective polynomial space. Hence, it is guaranteed that solutions of the 
optimization problems (^J, ([7J and ([8j exist. 

Remark 2.4. The optimization problem ^ has a well-known solution which can be found 
in |19[ Theorem 1.3.3]. The solutions for ([7]) and ([8]) given below can be considered as 
novel. 

Remark 2.5. The functional e(f) is by far not the only possible way to measure the 
space localization of a function. In the literature, there exist various other forms in 
this direction. In the Landau- Pollak-Slepian theory (cf. [7], |16| . |27j ) an optimization 
problem similar as in ^ is investigated. Results in terms of polynomial and exponential 
growth of polynomials can be found in the articles [6] and |14) . Results concerning the 
Shannon information entropy of orthogonal polynomials are summarized in the survey 
article [I]. 

In order to describe the optimal polynomials, we need the notion of associated and of 
scaled co-recursive associated polynomials. First of all, we know that the orthonormal 
polynomials p\ satisfy the following three-term recurrence relation (cf. f8j Section 1.3.2]) 

k +1 pi + i(x) = (x - ai)pi(x) - kpi(x), I = 0, 1, 2, 3, . . . (9) 

p-i(x) = 0, p (x) = —, 

with coefficients a/ € R and b\ > 0. 

Definition 2.6. For m £ N, the associated polynomials pi(x,m) on the interval [—1, 1] 
are defined by the shifted recurrence relation 

b m+l+1 pi +1 (x,m) = (x - a m+ i)pi(x,m) - b m+l p t -i(x , m) , I = 0, 1,2, . .. , (10) 
p_i(x,m) = 0, po(x,m) = 1. 

Further, for 7 G 1 and 5 > 0, we define the scaled co-recursive associated polynomials 
Pl(x, m, 7, 5) on [—1, 1] by the three-term recurrence relation 

b m+ i+iPi + i(x,m,-f,5) = (x - a m+ i)pi(x,m,7,6) - 6 m+ ^_i(x, m, 7, J), 

1 = 1,2,3,4..., (11) 
5x - a m - 7 



p (x, m, 7, 5) = 1, pi(x, m, 7, 5) 



m+l 
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The three-term recurrence relation of the co-recursive associated polynomials pi+i (x, m, 7, 5) 
corresponds to the three-term recurrence relation of the associated polynomials pi+i (x, m) 
except for the formula of the initial polynomial p\(x, m, 7, S). For c = 0, 7 = and 5=1, 
we have the identities pi(x, 0) = pi(x, 0,0, 1) = bopi(x). The polynomials pi(x,m) and 
Pl(x, m, 7, 5) can be described with help of the symmetric Jacobi matrix J™, < m < n, 
defined by 

/ a m b m+ i • • • \ 
b m +i «m+i b rn+ 2 • • • 






V 



b m +2 Om+2 





b n -2 «n-l b n -\ 







Jn-1 



a. 



(12) 



n ) 



If m = 0, we write J n instead of J^. Then, in view of the three-term recurrence formulas 



(10) and (11), the polynomials pi(x,m) and pi(x,m,j,5), I > 1, can be written as (cf. 
[E2 Theorem 2.2.4]) 

(13) 



and 



(x , m, 7, 5) = det I x 



Pi(x, m) = det{xli 
5 



1 



1-1 



'm+l-l)i 



T 771 

J m+Z-1 



7 
0,_i 



(14) 



where denotes the (I— l)-dimensional identity matrix and 0/_i the (I — l)-dimensional 
zero matrix. 

Next, we give a characterization of the functional e(P) in terms of the expansion coeffi- 
cients C[ of the polynomial P = X^/= n c lPl- 



Lemma 2.7. For the polynomial P(x) = cipi(x), we have 



1=0 



e(P) = c H 3 n c, 
e(P) = c H J™c, 
e(P) = c H J™c + (e(K) - a 



m ) I I 1 



if Peu n , 
if feC, 
if feC, 



with the coefficient vectors c = (cq, . . • , c n ) T and c = (c m , . . . , c n ) J . 



Proof. Using the three-term recurrence formula ^ and the orthonormality relation of 
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the polynomials pi, we get for P £ Tl n 

e(P) = / x \^cipi(x) w(x)dx = (y^cixpi(x)) [^cipi(x))w(x)d. 
1=0 J ~ 1 1=0 1=0 



/ + + b lPl-i(. x )j) (y~]cipi(xj)w(x)da 

1=0 1=0 
n n—1 

^2 a i\ c i\ 2 + ^{h+icici+i + k+lClCl+l) = c H J n c. 



1=0 1=0 

If Co = ... = c m _i = 0, we get the assertion for polynomials P in the space II™. If 
P G il^, then P has the representation 

(m— 1 \ n 

Pm{x) + ^2 e lPl( X ) + C lPl( X )> 
1=0 J l=m+l 

where the polynomial 1Z is given by lZ(x) = p m (x) + X/l=o e lPl( x )- Inserting this repre- 
sentation in the upper formula for e(P), yields the identity e(P) = (s(TZ) — a m )\c m \ 2 + 
c H J™c. □ 



Using the characterization of e(P) in Lemma 2.7 we proceed to the solution of the 
optimization problems (|6|, (J7| and ([8]). 

Theorem 2.8. The solutions of the optimization problems ^ and Q are given by 

n 

V n {x) = hi ^2pi(X n+ i)pi(x), (15) 

1=0 

n 

K{x) = K 2 Y,Pl-m{K-m+X,m)pi{x), (16) 
l=m 

Vn(x) = k 3 I K(x) + ^2 pi- m (\ 7 ^_ m+1 ,m^n^'R.)pi{x) J , (17) 

\ l=m+l / 

where pi(x,m) and pi(x, m, j-r, S-jz) denote the associated and the scaled co-recursive as- 
sociated polynomials as given in Definition 2.6 with the shift term 77^. := e(JZ) — a m and 
the scaling factor 5-ji := ||72.||^. 

The values A n +i, A™_ m+1 and A^_ m+1 denote the largest zero of the polynomials p n +\(x) , 
Pn—m+i{x,ni) and p n _ m +i 

(x,m,^n,5ii) in the interval [—1,1], respectively. The con- 
stants ki, K2 and K3 are chosen such that the optimal polynomials lie in the respective 
unit sphere and are uniquely determined up to multiplication with a complex scalar of 
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absolute value one. The maximal value of e in the respective polynomial space is given by 

M n := maxe(P) = X n+1 , 

Pes„ 

M™:= max. £ (P) = A™ m+1 , 
M? := max e(P) = A?L m , . 

Proof. We start out by determining the optimal solution V™ for the optimization prob- 
lem (J7|. The formula for the optimal polynomial V n follows as a special case if we 



set m = 0. First of all, Lemma 2.7 states that the mean value s(P) of a polyno- 
mial P(t) = X^ m QPz( x ) can be written as e(P) = c H J™c with the coefficient vec- 
tor c = (c m , • • • , c n ) T . Thus, maximizing e(P) with respect to a normed polyno- 
mial P £ is equivalent to maximize the quadratic functional c^J™c subject to 
l^| 2 = c m + c m+i + ' ' ' + c n = 1- ^ ^rT-m+i denotes the largest eigenvalue of the 
symmetric Jacobi matrix J™, we have 



c"J™c<A™ m+1 |c| 2 (18) 

and equality is attained for the eigenvectors corresponding to A™_ m+1 . Now, the largest 
eigenvalue of the Jacobi matrix J™ corresponds exactly with the largest zero of the 
associated polynomial p n _ m+ i(x, m) (cf. (U Theorem 1.31]). Using the recursion formula 



(10) of the associated polynomials pi(x,m) with c m = 1 the eigenvalue equation J™c 



A™_ m+1 c yields 

Q = pi- m (\™_ m+1 ,m), l = m,...n. 

Finally, we have to normalize the coefficients q, m < I < n, such that |c| 2 = 1. This 
is done by the absolute value of the constant K2- The uniqueness (up to a complex 
scalar with absolute value 1) of the optimal polynomial V n follows from the fact that 
the largest zero of p n _ m +i(x, m) is simple (see Theorem 5.3]). The formula for 



follows directly from the estimate in (jl S|) . 

states that in this case 



2.7 



We consider now the third polynomial space 11^. Lemma 
the mean value e(P) of P(x) = c m 1Z(x) + Yl?= m +i c iPi( x ) can ^ e written as e(P) = 
c^J™c + (s(TZ) — a m )|c m | 2 , with the coefficient vector c = (c m , • • • , c n ) T . Maximizing 
e(P) with respect to a polynomial P G is therefore equivalent to maximize the 
quadratic functional c H J™c + (e(TZ) — a m )\c m \ 2 subject to — l)|c m | 2 + |c| 2 = 1. 

Using a Lagrange multiplier A and differentiating the Lagrange function, we obtain the 
identity 

J™c + 7rc(c m ,0,--- ,0) T = X(5nc m ,Cm + i, - ■ ■ ,c n ) 
as a necessary condition for the maximum, where 77^. = e(TZ) — a m and 5-ji = ||7£||^. By 



the equation (14), this system of equations is related to the three-term recursion formula 



(11) of the scaled co-recursive associated polynomials pi(x, m,jiz, 5tz). In particular, 
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the value A corresponds to a root of p n - m+ \{x, m, 77^, 8n). Moreover, the maximum of 
c H J™c + 7rc|c m | 2 is attained for the largest root A = A^_ m+1 of p n ^ m+ i(x, m, jn, 8n) 
and the corresponding eigenvector 

c = K 3 (l,pi(X^_ m+1 ,m,jn,Sn), ■ ■ ■ ,Pn-m(A^_ m+1 , m,jn, <fo)) , 

where the constant K3 is chosen such that the condition (8n — l)|c m | 2 + |c| 2 = 1 is 
satisfied. The uniqueness of the polynomial (up to a complex scalar of absolute 
value one) follows from the simplicity of the largest root A^_ m+1 of the polynomials 
Pl{x, 171,711,811) (see [2j Theorem 5.3]). From the above argumentation it is also clear 
that the maximal value is precisely the largest eigenvalue A^_ m , 1 . □ 



3 Explicit expression for the optimally space localized polynomials 



Our next goal is to find explicit expressions for the optimal polynomials V n , V™ and 
derived in Theorem 2.8 To this end, we need a Christoffel-Darboux type formula for the 
associated polynomials pi(x, m) and pi(x, m, 7, 8). 

Lemma 3.1. Let pi(x,m) and m, 7, <5) be the associated and the scaled co-recursive 
associated polynomials as defined in (10) and (11). Then, for x 7^ y, the following 
Christoffel-Darboux type formulas hold: 

n 

^2pk(x)Pk-m(y,' 



, m 



k=m 



p n+ i(x)p n - m (y,m) - p n - m+ i(y,m)p n {x) p m -i(x) 

— On+l r o m . 

x — y x — y 

n 

^2pk{x)pk~ m (y,m,j,8) 

p n+ i(x)p„-m(y, m, 7, 8) - p n - m+1 (y, m, 7, 8)p n (x) 



k=m 



b n +i 

+ 



x-y 

Pm(x)((8 - l)y - 7) Pm-l(x) 



(19) 



(20) 



x-y 



x-y 



Proof We follow the lines of the proof of the original Christoffel-Darboux formula (see 
fj Theorem 4.5]). By (|9j) and ([ToJ) , we have for k > m the identities 

xp k (x)p k ^ m (y,m) 

= b k +ip k +i(x)p k -m(y, m) + a k pk(x)pk-m(y, m) + b k p k _ 1 (x)p k - m (y, m), 
ypk(x)p k - m (y,m) 

= b k+1 p k (x)p k - m+ i(y, m) + a k p k (x)p k ^ m (y, m) + b k Pk(x)pk-m-i{y, m). 
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Subtracting the second equation from the first, we get 

(x - y)Pk{x)Pk-m{y,rn) 

= h+i{Pk+i(x)Pk-m(y,m) - Pk(x)pk-m+i(y,™,)) 
- bk{pk{x)p k ^ m ^ 1 (y,m) - p k ^ 1 (x)pk- m (y,m)). 



Let 



Fk(x,y) = b k+1 



p k+1 (x)p k ^ m (y, m) - p k (x)p k - m+ i(y, m) 
x-y 



Then, the last equation can be rewritten as 

p k (x)p k -m(y,m) = F k (x,y) - F k -i(x,y), k > m, 



where F m -i(x,y) = —b m p m -i(x). Summing the latter from m to n, we obtain (19). 
Analogously, we get for the scaled co-recursive associated polynomials 

p k (x)p k - m (y, m, 7, 5) = G k (x, y) - G k -i(x, y), k>m + l, 
Pm{x)po{y,m,j,5) = p m (x), 



where 



G k (x,y) = b k+ 



Pk+i(x)pk-m(y, m, 7, 5) - p k (x)p k - m+1 (y, m, 7, 5) 



x-y 



G m (x, y) = b J^±M -Pmix){5y - a m - 7) fc > m + L 



x-y 



Then, summing from m to n, we get 

n n 

^2pk{x)pk-m(y,rn,j,5) = ^2 (Gk(x,y) - G k -i(x,y)) +p m (s 



k=m 



k=m+l 



b m+ ip m+ i(x) + p m (x)(5y - a m - 7) p m (x)(x - y) 

--G n [x,y) 1 

x — y x — y 

-G n (x, y) + PmixW-Dy-,) + ^p^x) 



x-y 



x-y 



Hence, we obtain formula (20). 



□ 



As a direct consequence of the Christoffel-Darboux type formulas in Lemma 3.1 we get 



the following explicit formulas for the optimal polynomials in Theorem 2.8 



Corollary 3.2. The optimal polynomials V n , V™ and Vl^ in Theorem 2.8 have the 
explicit form 



V n {x) = Kl&„+1 



p n+ l(x)p n (X n+1 ) 



x — X n 



+1 
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V n (x) = K 2 



T — \ m 



+ 



n—m+1 

Pm(x)((5 n - l)A^_ m+1 - 77?.) + 



t — \ n 

x A n-m+l 



where the constants k\, k 2) K3 and the roots A n +i A™_ m+1 and A^_ m+1 are given as in 
Theorem \2.8i. 



4 Optimally space localized polynomials for Jacobi expansions 

In this section, we will see that in the case of the Jacobi polynomials the generalized 
mean value e(f) is related to an uncertainty principle and that the optimal polynomials 
V n , V™ and V 1 ^ minimize the term vais(f) f° r the position variance of the uncertainty 
principle. 

In the following, the weight function w = w a p under consideration is the Jacobi weight 
function 

w aP {x) = (l-x) a (l+xf, x€[-l,l], *,P>-~- 

The corresponding orthonormal polynomials Pn (x) are called Jacobi polynomials and 
satisfy the differential equation L a gpn = —n(n + a + (3 + l)p^h'^\ where the second- 
order differential operator L a g is given as (cf. |28| Theorem 4.2.2.]) 

Lapf(x) = (1 - x 2 )^f(x) + (/3 - a + x(a + /3 + 2))^/(x). 

The next theorem states a well-known uncertainty principle for functions having an 
expansion in terms of Jacobi polynomials. 

Theorem 4.1. Let f £ C 2 ([-l, 1]) D L 2 ([-l, 1], w a p) such that \\f\\ Wa0 = 1. Further, let 

(a - p) + (a + + 2)e(f) ? 0. 
Then, the following uncertainty inequality holds: 

l~g(/) 2 It f f\ , (« + /3 + 2) 2 ( . 

T/ie constant ( Q +^+ 2 ) on f/ie rig/ii Zianc? side 0/ (21) is optimal. 
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Remark 4.2. Theorem 4.1 has been proven for ultraspherical expansions in |25j and was 
generalized to the Jacobi case in [18J. Hereby, the notation in [18] differs slightly from 



the notation above. For a more detailed discussion of Theorem 4.1 see also |3], [1], [9] 
and 1261. 



The terms 



var s (/) := 



2 > 



(22) 
(23) 



varp(/) := (-L aj3 f,f) Waf) 

in inequality (21) are called the position and the frequency variance of the function /, 
respectively. As the generalized mean value e(f), also the position variance varg(/) 
defines a measure for the localization of the function / at the boundary points of the 
interval [—1,1]. In particular, the more mass of the L 2 -density / is concentrated at 
the boundary points, the smaller the position variance varg(/) gets. The next Theorem 
shows that both measures are in principle equivalent. The only thing one has to take 
account of is that, in contrast to e(/), the position variance varg(/) does not differ 
between the two boundary points. Therefore, one has to restrict the set of admissible 
functions in the optimization problem: 



r 

>-'n 


:= {P G §n : 


e(P) > Ai}, 


rm 
L 'n 


:= {P G S- 


e(P) > Xi} 


c n 


:= {P G 


e(P) > Xx} 



where Ai = 2+a+p corresponds to the sole root of the Jacobi polynomial pf^\x) of 
degree 1. 

Theorem 4.3. If the sets C n , C™ and are nonempty, then 

arg min vais(P) = max e(P) = V n , 
Pec„ Pec n 



arg min var9(P) = arg max e(P) = V, 



rn 
n i 



arg min vais(P) = arg max e(P) = . 
p^ B Pe£* 

Hence, from all polynomials in the sets C n , C™ , C^, the optimal polynomials V n , V™ 
minimize the position variance varg. 

Proof. We consider the space variance varg as a function of A = s(f). We have 
var 5 (A) 



dvarg _ -2(A-Ai)A-2(l- A 2 ) _ -2(1 - AjA) 
~^A~ ( ' ~ (A- Ai) 3 ~ (A-Ai) 3 ' 



12 



5 Optimally space localized polynomials for Hermite expansions 



Therefore, the derivative 4r varg is strictly decreasing on the open interval (Ai,l) and 
strictly increasing on ( — 1, Ai). So, for P G C n , C™, maximizing s(P) yields the same 
result as minimizing the position variance vars(P). □ 

Remark 4.4. Whereas it can not be guaranteed that the sets £™ and C% are nonempty, 
the non-emptiness of the sets C n , n > 1, is a consequence of the interlacing property of 
the zeros of the Jacobi polynomials (cf. |28[ Theorem 3.3.2], |2j Theorem 5.3]). Namely, 
this interlacing property implies that e('P n ) = A n +i > A n > . . . > Ai. 



Remark 4.5. It can be shown that the uncertainty product v&rs(V n ) ■ varp^n) m (21) 
of the optimal polynomials V n is uniformly bounded by a constant independent of the 
degree n. Hence, the polynomials V n are not only well localized in space, but also in 
space and frequency. The quite technical proof can be found in [J]. 

Example 4.6. As a final example, we consider the orthonormal Chebyshev polynomials 
t n of first kind, i.e., the Jacobi polynomials Pn with a = (3 = — I and the weight 
function w a p{x) = 1. The orthonormal Chebyshev polynomials are explicitly given as 
(see [HI p. 28-29]) 

h( x ) = 7^F' t n (x) = ^cosnt, n > 1, 

where x = cost. The largest zero of the Chebyshev polynomials t n +i is given by A n +i = 
cos 2n+2 ( see 1^^' (6.3.5)]). The normalized associated polynomials t n (x,m), m > 1, 
correspond to the Chebyshev polynomials u n of the second kind given by (see [HI p. 28- 
29]) 

. , /2~sin(n+ l)t 

u n (x) = \ — , n > 0. 

V 7r smi 

The largest zero of the polynomials is given by A n +i = cos -^-^ ■ So, in the case of the 
Chebyshev polynomials of first kind, we get for the optimally space localized polynomials 
the formulas 

kir \ ki cosjn + l)t cos 
1 + 2 > cos cos kt = . 24) 

t[ 2n + 2 J * COSt -COS 

n ■ (k-m+l)n \ ( n-m+2 f \ t n+rn f \ 

Y Sm — +2 cosfct = 2^cos( -^t)co S (— t) ^ 

fc t^ S1I1 S=^+2 / n COst -™ST^[+2 

These optimal polynomials T n and T™{x) in combination with the Breitenberger uncer- 
tainty principle on the unit circle were intensively studied by Rauhut et al. in |23| and 




5 Optimally space localized polynomials for Hermite expansions 



As an example of an orthogonal expansion in a non-compact setting we consider the 
Hermite polynomials on the real line. The aim of this section is to construct polynomials 
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5 Optimally space localized polynomials for Hermite expansions 



having an expansion in terms of the Hermite polynomials that are optimally localized at 
the point x = 0. In the following, we will see that most of the theory of the previous 
sections can be applied also in this case, although with slight modifications. 

The Hilbert space under consideration is now L 2 (R, wh) with the weight function wh(x) = 
e~ x . The corresponding orthonormal polynomials (hi)f^Q, defining an orthonormal basis 
of L 2 (M, Wh), are called the (orthonormal) Hermite polynomials on R. 



As in Definition ( |2,1[ ), we introduce the polynomial spaces U n , n™ and n^, and the 
corresponding unit spheres § n , §™ and for the Hermite polynomials h n . The goal is, 
similar as in Section 2] to find those polynomials from S n , S™ and §^ that minimize the 



position variance 



,2 



var 5 (/) := / x'\f{x)\'e' x dx. (26) 



Since in the setting of the Hermite polynomials the calculations will be more complex, 
we will omit the case P G 1J^. Also, we will assume that n and m are even integers. The 
minimization problems then read as follows: 

Tin = arg min vars(P), (27) 

H™ = arg min var 5 (P). (28) 
Pes™ 



To get the explicit solutions, we will make use of the Laguerre polynomials p\ that 
form an orthonormal basis of the Hilbert space L 2 ([0, oo), w a ) with the weight function 

w a (x) = x a e~ x , a > —1. The orthonormal Laguerre polynomials pf^ and the Hermite 
polynomials hi are correlated by the following two formulas (see |13[ Section 4.6]): 

h 2 l(x) =p { - 1/2) (x 2 ), 1 = 0,1,2,..., (29) 

h 2 i+i(x) = xp\ 1/2 \x 2 ), 1 = 0,1,2,.... (30) 

The next Lemma gives a characterization of the position variance vars(P) in terms of 
the expansions coefficients of a polynomial P. 

n 

Lemma 5.1. For a polynomial P(x) = we get the formulas 

1=0 

var s (P) =cfj(-|)«c e + cfj(i) t _ lCoj if P G n n , 

m m 

var 5 (P)=cfj(-i)|c e + cfj(i)!_ 1 c , if Pen™, 

2 2 

with the coefficient vectors 

C e = (co, C2, . . . , C n ) T , C = (a, C3, . . . , C„_l) T , 

Ce = (Cm) C-m+2 ; ■ • • > Cn) , C Q = (c m -|-l , C m +3) • • • > C n — 1 ) , 
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and the matrices J(— j)| and J(^)|_ 1 corresponding to the Jacobi matrices of the as- 
sociated Laguerre polynomials pj ^\x,'y) an d Pi*\ x i**t)> respectively. 



Proof. Using the correlations (29) and (30) between the Hermite and Laguerre poly- 
nomials as well as the three-term recurrence formulas ([9j of the Laguerre polynomials 

ii n 
Pt 5 (a;) and p^ (x), we get for the polynomial P(x) = the formula 



1=0 



var s (P) 



1=0 



e x dx 



-1 



1=0 

n 
2 



C 2 l+1 X 3 p i l ~ 2 ) (X 2 ) ) ( °l k l ( X ) ) e ^ dx 
1=0 1=0 

5>(«k? 'i7W + oi - ^"* V) + ^ (_l) ptf V)) 



1=0 

1 



+ £ c^^p^s 2 ) + aJ^V) + ij^lfx 3 )) ) ( £? =0 cihi{x)\e~ x2 dx 



(§)_(§), 



-i 



+ ]T ca+i (6^^+3(1) + a; t; ^ +1 (x) + b^h^x))) (r:u c Mx))e- x2 dx. 

1=0 

Next, using the orthonormality relations of the Hermite polynomials hi, we can conclude 

for Pen„ 



.(I) 



(3), 



If c 



var s (P) =cf J(-|)nc e + cf J(|) ? _ lCo . 
c m _i = 0, we get the assertion for the polynomials P in n™. 



□ 



The solutions of the optimization problems (27) and (28) now read as follows. 



Theorem 5.2. The polynomials solving the minimization problems (27) and (28) are 
given by 



2 j 

K(s) = «i^P ; 2 (A| +1 ) /i 2i (x), 
1=0 

n 

_j_ / Is m 



(31) 
(32) 



15 
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where p l 2 (x,?j) denote the associated Laguerre polynomials with parameter a = —\. 

— (—-) 
The values A»+i an d A„ 2 _ m+1 denote the smallest zero of the polynomials p™+i {x) and 

Pn-m J^i y)> res P ec ti ve ly- The constants k% and K2 are chosen such that the optimal 
polynomials lie in the respective unit sphere and are uniquely determined up to multipli- 
cation with a complex scalar of absolute value one. The minimal value of vars(P) in the 
respective polynomial spaces is given by 

M n : = min vare(P) = As+i, 

Pe§„ 2 

m 

M™:= min var s (P) = A 2 r 



Proof. In the following, we will determine the optimal solution %™ for the minimiza- 
tion problem (28). The formula for the optimal polynomial % n follows as a special 
case when m = 0. By Lemma 5.1 the position variance vars(P) of a polynomial 

1 1 m m 

= Yl?=m c l^l( x ) can b e written as var s'(P) = c^J(— ^)| c e + c^ i{^)n_ 1 c with the 
coefficient vectors c e and c e given in Lemma 5.1 Hence, minimizing vars(P) with respect 
to a normed polynomial P 6 is equivalent to minimize the quadratic functional 

m m 

cf J(-i)| c e + c^J(i)|_ 1 c subject to |c e | 2 + |c | 2 = c 2 n + c 2 m+l + - ■ - + c 2 = 1. (33) 
The minimization problem ([5]) has a block matrix structure with the two matrices 

m m 

J(— i)| and J(i)n_ v To solve the problem, we have to determine the smallest eigen- 

2 2 

m m 

value of both matrices. The eigenvalues of the two matrices J( — 1)| and J(i)ni corre- 

^2 ^ 2 1 

spond exactly with the zeros of the Laguerre polynomials p n _ 2 m (x, and p„ 2 _ m (x, ™ ). 

~~ 2 1-1 ~~2 ~ 

By a functional analytic method based on the Hellmann-Feynman Theorem (see |12j , |13[ 
Section 7.4] and |5]) it follows that the smallest zero of the associated Laguerre polyno- 
mials p^ a \x,m) is an increasing function of the parameter a. This result together with 

(_!) 

the fact that the smallest eigenvalue of p n 2 m (x, y) is strictly smaller than the one 
°f Pn-m { x i y) ( see ^ ne interlacing property of the orthogonal polynomials [28,, Theorem 

2 

m 

3.3.2]) implies that the matrix J(— i)| is the one with the smallest eigenvalue. 



2J 2 

V L\ 11 V ' IV :f U IA OXA-LClllv^O I ■ im I1W II ' L\ V '.L 1 IV -MIV VI '1 1 1 i< V 11 .J.,\ *f I r, 



Therefore, if A _ L1 denotes the smallest eigenvalue of the Jacobi matrix J(— ^)| , we 
get 

1 m 1 77i m 

cf J(-^)Ic e + cf J(-)I_ 1 c > A J (|c e | 2 + |c | 2 ) (34) 

2 2 2 2 1 —5— + 1 

m 

and equality is attained for the eigenvectors corresponding to A 2 _ m . Using the recur- 
I I 

sion formula (10) of the associated polynomials p l 2 (x,^j) with c m = 1 the eigenvalue 
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equation J( — I) J c 



A 2 

n — m I -i 
2 ~ rL 



c e yields 



(-21 
C21+1 



0, 



2 

n — m 
2 



m \ 
+1' 2 )> 



m n 
2 ' • • • ' 2 ' 



m n 
2 ' • • • ' 2 



1. 



Finally, we have to normalize the coefficients q, m < I < n, such that |c e | 2 + |c G | 2 = 1. 
This is done by the absolute value of the constant K2- The uniqueness (up to a complex 
scalar with absolute value 1) of the optimal polynomial V n follows from the fact that 

i( x > If) * s s i m pl e ( see El Theorem 5.3]). The 

□ 



the smallest zero of the polynomial p 



n — in 
2 



formula for M™ follows directly from the estimate in ( 34 ) . 



Using once again the relation (29) between the even Hermite polynomials and the La- 



guerre polynomials and the Christoffel-Darboux type formulas ot Lemma |3.1| we get the 



following explicit formulas for the optimal polynomials in Theorem 5.2 



Corollary 5.3. The optimal polynomials 7~L n and 7-L™ in Theorem 5.2 have the following 
explicit form: 

n, t \ J- D (-Dm a h n+ i(x) 



Aa+i 



K2- 



\ 2 

'* n — m 



+ 1 



where the constants K\, k<i and the roots Aa_|_i, 



are given as in Theorem 



5.2 



6 Construction of polynomial filters for the detection of peaks in periodic signals 



In this final section, we will give some examples on how the optimal polynomials of 
Section [4] can be applied as filters for the detection of peaks. To this end, we consider 
continuous 27r-periodic signal functions / E L 2 ([—tt,tt)) and trigonometric polynomial 
filters h G 

Our goal is to find trigonometric polynomials h which are well suited to work out the 
peaks of the signal /. Since the filtering operator defined in ([TJ acts as a convolution 
operator on /, most of the mass of the polynomial h has to be concentrated at the 
point t = in order to filter out the peaks of the signal /. If we further assume that h 
is even, i.e., h{t) = h(—t), then /i(arccos(x)) is defined on the interval [—1,1] and is a 
polynomial of degree n in the variable x = cost. Moreover, if the polynomial /i(arccos(x)) 
is localized at x = 1, then the trigonometric polynomial h is localized at t = 0. Therefore, 
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6 Construction of polynomial Biters for the detection of peaks in periodic signals 



the optimally space localized polynomials 7-^ (cost) of Theorem 2.8 are natural choices 
for polynomial filters in peak analysis. 



Experimenting with different weight functions w in Theorem 2.8 and Corollary 3.2, it 
is possible to construct a whole bunch of well-localized polynomial filters with different 
properties. We give here just some easy examples using the Jacobi weight function w Q p. 

i, j3 = — i, we get the two filter kernels 

cos(n + l)i 



For a = — i, /3 



and a 



/#)(*) :=T n (cost) = CW 



cos t — COS 

1 

COS 



2n+2 



(35) 
(36) 



where the constants Cn and C n denote normalizing factors such that the respective 

polynomials are norm ed i n the L 2 -norm. The optimal polynomial /in (i) was already 

(i) 



4.6 



computed in Example 
are known as Rogosinski 
known as cosine windows (cf. [llj). 



In approximation theory, the trigonometric polynomials h n 
ternels (cf. |17| p. 112-114]), in signal analysis they are well- 



The Rogosinski filter h t 



(i) 



The polynomial filter h, 



(2) 



Figure 1: The polynomial filters and of degree 6. 

The filter h^f 1 (t) is computed in the same way as hn (t) using the explicit representation 
of the Chebyshev polynomials of third kind (for the definition, see [HI Section 1.5.1]). 
Figure [l illustrates that compared to the filter the trigonometric polynomial h 
has a wider peak at t = but less mass at the ends t = tt and t = —tt. This is due to 

(2) 

the fact that in the case of the filter h n we optimize over all polynomials P £ U n with 



j\p[x)\\l 



i . _ i 
x) 2 (1 + x) tax 



i.e. the particular optimization problem favours polynomials that have more mass con- 
centrated at x = 1. 
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2 4 0.6 0.8 1 

Optimal polynomial, n = 12 




0.2 D 4 0.6 0.8 1 

Optimal polynomial . n=32 




0.4 0.6 
2. Filtered signal 



Figure 2: Filtering a noisy signal with the optimal polynomial filter hffl ■ 



If the peaks of the signal / lie on a low-frequency carrier signal or if some baseline 
correction has to be done, it is reasonable to additionally filter out the low frequencies 
of /. In this case, the polynomial filter h has to be restricted to a frequency band 



[m, n] C N. The corresponding optimal filters are given in equation (16). In case of the 
Chebyshev polynomials U(x) of first kind, the optimal band-limited polynomials are 



given in Example 4.6 



as 



cos( 



n— m+2 



t) cos(™i) 



cos t — cos 



(37) 



n—m+2 



Finally, we consider polynomial window functions h with the additional property h(ir) = 
h{—Ti) =0. If h is assumed to be a symmetric trigonometric polynomial, then h has the 
form hit) = g{t)(l + cost), where g is a symmetric trigonometric polynomial of degree 
n — 1. Now, as above, we can play the same game with the polynomial g, i.e., we can 
insert optimally space localized polynomials as candidates for g. Using Jacobi weights 
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with a = — |, (3 = | and a = |, /3 = |, we can introduce the filter kernels 

cos t — X n 

h^t):=0^ S ^±m '+"". , (39) 
sin £ cos £ — cos ^h- 

where A n denotes the largest zero of the Jacobi polynomial p n 2 ' 2 . The filter h-h is 
defined such that the functional x\P(x)\ 2 (l + x) 2 dx is maximized over all polynomials 

P under the constraint j\ \P(x)\ 2 (l+x) 2 dx = 1. The filter polynomial hffl is well-known 
in signal analysis under the name Hann window (see In fact, computing the Fourier 

coefficients of hj% , one gets hn (j) = c cos 2 ( 2 n+2 ) ^ or J = — n > • • • J n - 
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